Triaxial super deformation in 40 Ar 
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Abstract 

Superdeformed (SD) states in 40 Ar have been studied using the deformed-basis antisymmetrized 
molecular dynamics. Low energy states were calculated by the parity and angular momentum pro- 
jection (AMP) and the generator coordinate method (GCM). Basis wave functions were obtained 
by the energy variation with a constraint on the quadrupole deformation parameter /3, while other 
quantities such as triaxiality 7 were optimized by the energy variation. By the GCM calculation, 
an SD band was obtained just above the ground state (GS) band. The SD band involves a K n = 2 + 
side band due to the triaxiality. The calculated electric quadrupole transition strengths of the SD 
band reproduce the experimental values appropriately. Triaxiality is significant for understanding 
lowdying states. 

PACS numbers: 
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Numerous superdeformed (SD) bands have been identified and studied in various domains 
of the nuclear chart. An example of recent progress in this research area is the discovery of 
the SD states in a very light mass region (A ~ 40) such as 36 Ar and 44 Ti [3J. A 

striking feature of these experimental works is that the SD band members were assigned from 
band heads up to high spin states. Discrete 7 transitions linked those SD bands to spherical 
or normal-deformed states, which established the excitation energies and the spin-parities 
of the SD bands. 

Theoretical microscopic studies of deformed states in the A ~ 40 region have been per- 



formed based on the shell model j^, 5], mean-field models 



-|8fl, and the antisymmetrized 



molecular dynamics (AMD) 9l-|l4|. Calculations with no assumption for the axial symme- 
try suggest the systematic triaxial deformation of the SD states. For example, the AMD 
calculations suggest that SD bands in 40 Ca [n]] and 44 Ti {9] are constructed from triaxial 
shapes and their side bands exist. The triaxiality of superdeformations in this mass region 
has also been discussed in the three-dimensional coordinate- mesh Hartree-Fock (HF) 6|. 
Triaxiality is also an important feature of normal-deformed states in the A ~ 40 region. 



Established triaxia 



states are in 40 Ca, in which an observed K 11 



2 + band with low 



normal-deformation 



excitation energies 2 1 is interpreted as a side band of a K w = + band caused by the triaxial 
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161 ] . Therefore, triaxiality is one of the key properties to clarify 



the structure of normal- deformations and superdeformations in A ~ 40 nuclei. However, 
theoretical studies taking into account triaxiality are still limited due to the higher numerical 
costs than axially symmetric calculations. 

Quite recently, an SD K* = + band built on the J 7 * = 0+ state (2.12 MeV) in 40 Ar was 
experimentally identified up to the J 71 " = (12 + ) state by the 26 Mg( 18 0, 2p2n) 40 Ar reactions 



171 ]. which were the lightest SD states ever in N 7^ Z nuclei. Superdeformations in N 7^ Z 



nuclei in this mass region are frontier topics. 

To discuss triaxiality, features of unnatural spin-parity [71 = (— ) J+1 ] states are helpful 
data. In 40 Ar, a state at 4.23 MeV was assigned to an unnatural parity J n = 3 + state by the 
40 Ar[a, a') reaction 18[. This J n = 3 + state decays to the J n = 2 + state in the SD band 
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20|. The .r = 2+ state at 3.92 MeV just below the J 71 " = 3 + at 4.23 MeV also decays to 



ihe SD band member J w = + , and the E2 transition strength takes a non- negligible value 
2o| . These observations show that the states J n = 2 + and 3 + are candidates for members 
of the side band in the SD band, which can be discussed in relation to the triaxiality of the 
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superdeformation. Nevertheless, neither theoretical studies of the SD state nor the triaxiality 
of 40 Ar have progressed. 

The purpose of this study is to investigate the structure of the SD band in 40 Ar focusing on 
triaxiality. Competitions of deformations favored by protons and neutrons are also discussed. 

To obtain the wave functions in low-lying states, the parity and angular momentum pro- 
jection (AMP) and the generator coordinate method (GCM) with deformed-basis AMD wave 
functions were performed. A deformed-basis AMD wave function |<3>) is a Slater determinant 
of triaxially deformed Gaussian wave packets, such that 

|$) = A\<Pi,<p 2 ,---,<Pa), (1) 

\<Pi) = |&)®lx<)®fa), (2) 
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\Xi) = Xi\ t> + Xi 1 I), ki> = k) or k>, ( 4 ) 

where A is the antisymmetrization operator and \ipi) are single-particle wave functions. The 
variables \xi), and \ri) are a spatial part, a spin part, and an isospin part, respectively, of 
every single-particle wave function The values u = (u x , u y , u z ) are real parameters that 
indicate the width of the Gaussian single-particle wave function, and Zj = (Z ix , Z iy , Z iz ) are 
complex parameters that indicate centroids of the single-particle wave functions. The width 
parameters v are common to all nucleons. The complex parameters x\ an d Xi represent 
spin directions. Basis wave functions of the GCM were obtained by the energy variation 
with a constraint on the quadrupole deformation parameter of a total system after the 
projection to positive parity states. Variation parameters were u, Z i; and x\ ■ The isospin 
part of each single-particle wave function was fixed as a proton or a neutron. The Gogny 



D1S 
Refs. 



brce was used as an effective interaction. Details of the framework are provided in 
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Calculated values of energy, triaxiality 7, and harmonic oscillator (H.O.) quanta AA^h.o. 
for protons and neutrons measured from the lowest allowed states are shown as functions of 
the quadrupole deformation parameter (3 in the upper, middle, and lower panels, respectively, 
of Fig. [TJ As for the energy curves, energies of the positive parity states and J n = + and 
2 + states, which were projected from each basis wave function, are plotted, i^-mixing was 
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FIG. 1: (upper) Energies projected to positive parity states, J w = + and 2 + states, (middle) the 
triaxiality 7 for protons and neutrons, and (lower) the H.O. values relative to lowest allowed states 
are plotted as functions of the quadrupole deformation parameter j3 of the total system. 

performed for the J 71 " = 2 + states. Those energy curves have two local minima at /3 ~ 0.2 
and 0.5, which are labeled GS and SD minima, respectively, because the GCM results show 
that these minima correspond to the GS and the SD states, respectively, as explained later. 
In the energy curves projected to positive parity states, the energy gap between the GS and 
SD minima is approximately 6.7 MeV, while the gap decreases to 2.5 MeV after the AMP 
to ,r = 0+ states. 

The values of triaxiality 7 for protons and those for neutrons are similar in the entire 
f3 region (the middle panel of Fig. [T|). In the small /3 region, the system forms an oblate 
shape with 7 ~ 60°, while in the large (3 region around the SD minimum, it forms a triaxial 
deformation with 7 ~ 10°. Density distributions of the obtained wave function at the SD 
minimum [(^,7) = (0.478,11.2°)] are shown in Fig. [2j Because of triaxial deformations 
around the SD minimum, each wave function contains different K components, and two 
j7T _ 2+ states are obtained after diagonalizing those different K components, i.e., K- mixing 
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FIG. 2: Density distributions of the obtained wave functions with (/5, 7) = (0.478, 11.2°) projected 
onto the xy, zy, and zx planes. 

(the upper panel of Fig. 

The relative H.O. quanta AiVg and AN^ Q for protons and neutrons, respectively, as 
functions of the quadrupole deformation parameter (3 (the lower panel of Fig. [I]) show the 
structural changes along (3 in terms of particle-hole excitations across major shell gaps. The 
AiV^'o are defined as 



AN n ' u 

ZAiv H.O. 

= E ( E £ + 



3 -A 

2 



-iV" 



(5) 



Here m is the nucleon mass, and the frequency uj = (u x ,u y ,uj z ) of the H.O. was chosen to 
be the deformed one uj = —u with the width parameters v of each wave function. A^ in and 
N^ in are H.O. quanta of the lowest allowed state for protons and neutrons, respectively, and 
in the case of 40 Ar, (iVj£ in , A^ in ) = (26, 36). In the small (3 region, the relative H.O. quanta 
AA^h.o. take an approximately zero value and indicate no particle-hole excitation for both 
protons and neutrons. On the other hand, in the f3 ^ 0.4 region, the values of AiVH.o. jump 
to two for both protons and neutrons. This shows that the SD minimum has 2hu excited 
configurations for both protons and neutrons, which is a total of 4hu excited configurations. 
In (3 ~ 0.3, ( ATVjj , AiVg ) = (1, 1) components are obtained. 

To investigate the single-particle behavior in detail, energies and parities of single-particle 
orbits of the obtained basis wave functions were calculated by transforming the deformed- 



basis AMD single-particle wave functions to the HF single-particle orbits 
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21] . Since most 
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FIG. 3: Single-particle energies for protons (upper) and neutrons (lower) are plotted as functions 
of the quadrupole deformation parameter of the total system. 

of the obtained wave functions have approximately parity symmetric shapes, extracted 
single-particle orbits are almost pure-parity orbits. Figure E] shows single-particle ener- 
gies of positive- and negative-parity orbits for protons and neutrons. Two orbits degenerate 
because of the time reversal symmetry except for ~ and 0.3. In the small region with 
< 0.25, the single-particle orbits are normal order for both protons and neutrons, which 
shows that configurations around the GS minimum are roughly [(s<i)~ 2 ] 7r [(p/) 2 ]^ relative to 
the ground state of 40 Ca. With an increase of deformation 0, the highest positive-parity 
orbits in the s<i-shell disappear at ~ 0.3 and negative-parity orbits come down for both 
protons and neutrons, which corresponds to the orbit inversion when the ^/-orbits intrude 
into the s<i-orbits in the large region. Consequently, configurations in the large region 
with > 0.35 around the SD minimum are roughly [(sd)~ 4 (pf) 2 ] n [{sd)~ 2 {j)f) A ] v . 

Figure H] shows the experimental and theoretical excitation energies obtained by the 
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FIG. 4: Level scheme of 40 Ar is shown. In experimental data, the GS, the SD band (SDO), 
candidate states of the side band of the SD band (SD2), and other low-lying J 71 " < 4+ states are 



listed. The experimental data are taken from Refs. 



17 and 



24. 



GCM. The K n = 0+ GS band, the K* = 0+ SD band (SDO), and the K* = 2+ SD band 
(SD2) are obtained by the GCM calculation. The energy levels of the GS and SDO bands 
agree with the experimental data. The observed J n = 2 + (3.92 MeV), 3 + (4.23 MeV), and 
(2,3,4)+ (5.17 MeV) states are tentatively assigned as members of the SD2 band in this 
study because their excitation energies are similar to calculated ones of the SD2 band and 
the candidates states decay to the members of the SDO band. Dominant components of 
the GS band are the deformed-basis AMD wave functions around the GS minimum on the 
energy curve, and those of the SDO and SD2 bands are the basis wave functions around the 
SD minimum with the triaxial deformations. Owing to the triaxial deformations around the 
SD local minimum, the SD2 band appears. The members of the SDO and SD2 bands are 
dominated by the \K\ = and 2 components, respectively, projected from a single deformed- 
basis AMD wave function at (/?, 7) = (0.478,11.2°) up to high spin states. Both the SDO 
and SD2 bands show rigid rotor-like behavior. The calculated level spacings in the SDO 
band members show the rotational band spectra and match experimental data, and those 
of the SD2 band members are consistent with experimental data. The energy gap between 
the J n = + D0 and the J n = 2 + D2 states is slightly overestimated, and it may indicate that 
the triaxiality of the SD states is underestimated. 
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The level spacings of the J n < 6 + states in the calculated GS band are rotational, but the 
level spacings between the J 71 = 6 + and 8 + are smaller than those in the case of a rotational 
band (Fig. SJ). The 2hu [(AN£ Q , AN&_ Q ) = (1,1)] components mix in the J 71 " = 8+ state; 
that breaks the rotational spectra. The calculated GS band deviates from the experimental 
one and is more rotational. The experimental GS band may be reproduced by adding 
sufficient basis wave functions such as (AiVg , AiVg ) = (1,1) components. The observed 
J w = 2 + (3.21 MeV) and 2 + (3.51 MeV) states were not obtained perhaps because basis 
wave functions with a total of 2hu configurations are insufficient in the present basis set. 

Table [I] shows theoretical and experimental values of electric quadrupole transition 
strengths B(E2) in the Weisskopf unit. The B(E2) values for the intraband transitions 
of the SDO and SD2 bands take larger values than those of the GS band, which reflects 
the large deformations of the SDO and SD2 bands. The B(E2) values for the intraband 
transitions of the GS and SDO bands are consistent with experimental data. Since the SDO 
and SD2 bands are approximately regarded as the K w = + and 2 + bands associated with 
the triaxial superdeformation, the interband transitions between them have non-negligible 
values. The observed J w = 2 + (3.92 MeV), 3+ (4.23 MeV), and (2,3,4)+ (5.17 MeV) de- 
cay to the SDO band members. The agreement of the theoretical and experimental B(E2) 
values of the 2 + D2 — > + D0 transition supports the band assignment of the SD2 band for 
the J n = 2j (3.92 MeV) state. B{E2) ratios of B(E2; 2 + D2 0^ D0 )/B(E2; 2 + D0 + D0 ) 
and B{E2\ 2 + D2 — > 2 + nn ) / B(E2; 2 + T-, n — > + D0 ) are consistent in the case of 7 ~ 10° with 



the Davydov-Filippov model 



25j, which is a schematic triaxial rotor model. To confirm 



the proposed band assignment of the SD2 band and triaxial deformation of the SD states, 
more experimental data for intraband transitions of the SD2 band and interband transitions 
between the SDO and the SD2 are required. 

To discuss shell effects of protons and neutrons in the SD states of 40 Ar, triaxial aspects 
of the superdeformation in 40 Ar[(Z, N) = (18, 22)] are compared with those in the N = Z 
isotope ( 36 Ar[(Z,iV) = (18,18)]) and isotone ( 44 Ti[(Z,iV) = (22, 22)}Y The deformed- 



13| and 44 Ti 



basis AMD and HF calculations suggest that the SD states in 36 Ar 
form prolate and significantly triaxial shapes, respectively. The theoretically suggested 
triaxiality 7 values of the SD states in 44 Ti are approximately 25° Q] by the deformed- 
basis AMD, which is consistent with the HF 6|. The particle-hole configurations of the SD 
states in 36 Ar and 44 Ti are considered to be [(sd)~ 4 (p/) 2 ] and [(sd)~ 2 (pf) 4 }, respectively, for 
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TABLE I: Electric quadrupole transition strengths are shown in the Weisskopf unit. Jj and Jf 
indicate initial and final states, respectively. The blanks under theoretical values indicate that 
calculated B(E2) values are almost zero. As for the transitions labeled " — " under experimental 



values, 7 decays have been observed, but mu 



The experimental data are taken from Ref. 
as the 4g D2 . 
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tipolarities and strengths have not been measured. 
The observed (2,3,4)+ state at 5.17 MeV is noted 
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both protons and neutrons, which indicates that the former configuration favors a prolate 
shape and the latter one induces a triaxial shape. The SD states in 40 Ar have dominantly 
the [(sci)~ 4 (p/) 2 ]^[(s(i)~ 2 (p/) 4 ] ;/ configuration, which shows an admixture of shell effects of 
these configurations. The triaxiality 7 values of protons and neutrons of the SD states 
in 40 Ar are approximately 10° (Fig. UJ), which is the intermediate value between 7 ~ 0° 
and 25° for the SD states in 36 Ar and 44 Ti, respectively. Protons and neutrons may favor 
similar deformations to increase overlap and gain attractive forces between them. Thus, the 
triaxiality of the SD states in 40 Ar is understood by a competition of a prolate shape favored 
by proton shell effects for Z = 18 and a significantly triaxial shape induced by neutron shell 
effects for N = 22. 

In conclusion, the structure of the SD states in 40 Ar has been investigated with the 
deformed-basis AMD and the GCM. The SD states in 40 Ar form triaxial shapes. Owing 
to the triaxial deformation, the K n = 2 + SD band exists. Calculated B(E2) values for 
intraband transition of the GS and SDO bands are consistent with experimental data. The 
triaxiality 7 of the SD states is understood by the competition of deformations favored 
by their proton and neutron structures. To understand details of nuclear structures, it is 
necessary to take into account triaxiality. 

Numerical calculations were conducted on the RIKEN Cluster of Clusters (RICC) and 
the supercomputer SX9 in the Research Center for Nuclear Physics, Osaka University. The 
authors thank Dr. N. Itagaki in the Department of Physics, Faculty of Science, The Univer- 
sity of Tokyo and Dr. T. Yoshida in the Center for Nuclear Study, The University of Tokyo, 
for fruitful discussions. 
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